ad  6059511 


REG.  NO _ 

Loa 


COKVECTION  HEAT  TRANSFER  CCKFFICIENTS 
IN  DE  LAVAL  NOZZLES 


TECHNICAL  REPORT  NO.  HS-3 


C' 


COW  OFTT 


hard  copy 
microfiche 


$. 

$. 


» 


mY  CONTRACT  N6-onr-25l  TASK  ORDER  6 
(NR-090-10if) 


Department  of  Mechanical  Engineering 
Stanford  University 
Stanford,  California 


October  1,  1951 


Prepared  Byi 
R.  K.  Pefley 


Approved  byS 
A,  L«  London 
Project  Supervisor 


ABSTRACT 


A/ 

The  objectives  of  this  program  are  to  determine  the  local 
convection  heat  transfer  coefficients  in  the  convergent  portion 
and  throat  sections  of  de  Laval  nozzles  2  and  to  determine  the  ef¬ 
fect  on  the  local  heat  transfer  coefficients.,  particularly  at  the 
throat,  of  varying  the  contour  in  the  converging  section  of  the 
nozzle. 

Three  two-dimensional  nozzles  -^^iasprlbod  herCT-ftei*  'as-ffTy 
th  ■  inad  i  w  TUyi  f  In,.  and..,i»aBn1 

-in£jgi.Q55j?fl€ot?l‘onal^areas^aa^dea'C.ilbgd  ±n  PlgT'lbN,,  have  been  tested , 
Heat  transfer  coefficients  have  been  obtained  for  both  the  flat  and 
contoured  walls., 

'  /T" 

The  results  do  not  correlate  when  using  a  length  dimension 
based  on  hydraulic  diameter,  as  conventionally  employed  for  estab¬ 
lished  pipe  flow. 

The  results,  for  both  the  flat  and  contoured  walls,  do  corre¬ 
late,  when  using  a  length  dimension  measured  from  the  nozzle  entrance 
along  the  walls ,  However ,  the  heat  transfer  coefficients  are  higher 
than  for  turbulent  boundary  layer  flat -plate  flow.  This  appears 
to  be  partially  due  to  the  technique  used  in  obtaining  the  appar¬ 
ently  local  coefficients,  and  partially  due  to  the  negative  pressure 
gradient  effect, 

A  significant  result  of  these  tests  is  that  the  throat  coeffi¬ 
cients,  for  the  same  mass  flow  rates  through  all  the  nozzles  III, 

IV,  and  V,  agree  within  t  15%,  indicating  that  the  markedly  dif¬ 
ferent  throat  approach  contour  variations  shown  in  Fig,  2  do  not 
produce  significant  variations  in  the  nozzle  throat  coefficients, 

A  theoretical  solution,  predicting  the  local  heat  transfer  co¬ 
efficients  along  the  flat  wall  of  nozzle  III,  has  been  made  l^y  ^- 
troducing  the  pressure  gradient  of  nozzle  III  into  the  von  Karman 
boundary  layer  momentum  equation.  The  results  of  this  work  are  com¬ 
pared  to  the  experimental  findings,  and  are  utilized  to  interpret 
the  experimental  results , 
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NOMENCLATURE 


English  Letter  Syfabols 

2 

A  -  heat  transfer  area,  ft 
a  -  length, ft 

C  -  specific  heat  of  solid  materials,  Btu/.'(lb  ®P) 
c  -  specific  heat  of  fluid  at  constant  pressure,  Btu/(3)h  ®F) 

D  -  hydraulic  diameter  of  nozzles,  r^^,  ft 
G  -  mass  flow  velocity,  lbs/!(hr'"ft^) 

h  -  unit  convective  heat  transfer  coefficient,  Btu/(hr  '“J*  ft^) 
k  -  thermal  conductivity,  Btu/(hf  ft^  ®F./f t) 

L  -  total  nozzle  length  to  throat,  ft 
m  -  arbitrary  length  dimension,  ft 
p  -  pressure,  Ibs/ft^ 
q  -  heat  flow  rate,  Btu/hr 
r  -  radius  from  center  of  plug,  ft 
t  -  temperature,  °F 

q 

V  -  total  volunle  of  copper  plug,  ft 

X  -  distance  from  nozzle  entrance  measured  along  the  wall 
profile,  ft 

Greek  Letter  Symbols 

cA  -  thermal  diffuslvity,  ft^/hr 
S  -  boundary  layer  thickness,  ft 
^  -  finite  differenthci 
P  -  density,  Ibs/ft^ 


ill. 


0  -  time,  hr 

-  viscosity  lbs /(hr  ft  ) 

Dimensionless  Groups 

St  -  Stanton's  Number  h/G  cl 

P 

Pr  -  Prandtl's  Number  c^/k 

Nu  -  Nusselu's  Number  h  x  /k  or  hD/li: 

Re  -  Reynold's  Number  G  x^ic^  or 
M  -  Mach  Number 
Subscripts 

b  -  denotes  bakelite  matrix 

aw  -  denotes  adiabatic  wall  temperature 

0  -  denotes  stream  stagnation  conditions 
0  -  denotes  thermal  plug  variable  with  respect  to  time 

1  -  denotes  conditions  between  thermal  plug  and  gas  stream 

2  -  denotes  conditions  between  thermal  plug  and  bakelite  matrix 

T  -  denotes  total 

r  -  denotes  any  constant  radius 
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INTRODUCTION 


A  detailed  description  of  the  purpose,  proposed  test  program, 
test  apparatus  and  some  test  results  has  been  made  in  two  previous 
reports,  references  (d)  and  (e) ,  To  facilitate  reading  this  report 
these  topics  will  be  briefly  summarized o 


Ob.iectives  of  Program 

To  determine  the  local  convection  heat  transfer  coefficients 
in  the  convergent  section  of  a  series  of  de  Laval  nozzles , 

To  discover  whether  or  not  modification  of  the  curvature  of 
the  nozzle  wall  approaching  the  throat  section  will  have  significant 
influence  on  the  convection  coefficients  at  the  throat « 

To  formulate  correlations  of  results  that  may  be  useful  to 
rocket  designers  and  others  concerned  with  heat  transfer  in  nozzles. 


Test  Program 

To  test  a  series  of  five  two  dimensional  nozzles,  the  only 
variables  being  the  throat  approach  contour  and  th'  mass  flow  rate. 

To  make  further  tests  and  modifications  to  the  test  apparatus 
as  dictated  by  the  results  obtained  from  these  experiements , 


Test  Apparatus 

The  series  of  nozzles  are  made  from  bakelite  to  minimize  the 
undesirable  heat  loss  and  yet  have  sufficient  strength  to  prevent 
distortion  or  leakage.  They  are  two  dimensional,*  i.e,,  one  pair 
of  opposite  walls  is  flat  and  one  pair  is  curved. 

Copper  plugs,  3/8  inches  in  diameter,  are  buried  in  one  flat 
wall  of  the  nozzle  and  one  curved  wall,  and  are  spaced  uniformly 
along  the  axis  of  the  nozzles ,  Their  faces  are  machined  flush 
with  the  faces  of  the  bakelite  (See  Fig,  2), 

Local  heat  transfer  coefficients  are  obtained  by  a  t.  isient 
technique,  •'■he  copper  plugs  are  heated  by  a  hot  iron  pressed  agains 
their  back  side  —  then  cooled  by  the  high  velocity  gas  stream 
sweeping  their  front  side. 


Previously  Reported  Results 


Local  heat  transfer  coefficients  were  reported  for  the  contour 
walls  of  nozzles  II  and  III .  Sufficient  data  was  not  available 
to  permit  correlation  in  any  suitable  form. 

The  limited  results  tentatively  indicated  that  a  considerable 
variation  in  the  local  heat  transfer  coefficient  at  the  throat 
might  be  experienced  by  modifying  the  throat  approach  contour. 


Purpose  of  Current  Report 

(a)  To  analyze  results  obtained  in  terms  of  available 
theory  and  conventional  correlations , 

(b)  To  further  consider,  in  the  light  of  additional  data, 
the  tentative  conclusion  that  contour  variation  in 

the  convergent  section  would  have  a  considerable  effect 
on  the  throat  heat  transfer  coefficient . 

(c)  To  describe  a  theoretical  attempt  to  predict  the  local 
coefficients  in  a  nozzle,  neglecting  the  effect  of  wall 
curvature  but  considering  the  effect  of  pressure 
gradient , 
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SUMMARY  OF  RESULTS  AND  CONCLUSIONS 


For  this  method  of  testing  and  v^ith  the  nozzle  contour  var¬ 
iations  described  in  the  body  of  the  report,  the  variation  in 
throat  coefficient  as  a  function  of  modification  of  wall  curva¬ 
ture  is  not  great  enough  to  single  out  any  one  contour  as  being 
preferred.  The  variation  in  contour  produces  considerable  local 
variation  in  coefficient  due  to  wall  curvature,  but  this  effect, 
along  with  the  effect  of  pressure  gradient  variation,  does  not 
persist  into  the  throat  section. 

The  data  obtained  clearly  demonstrates  that  prediction  of 
heat  transfer  coefficient*’  in  the  convergent  portion  of  a 
de  Laval  nozzle  should  be  made  on  a  flat  plate  or  pipe  entrance 
type  of  correlation,  (Contrast  Figs,  Va,  b  and  5a,  b) , 

If  the  nozzle  is  long  (L/D  large  based  on  throat  hydraulic 
diameter)  it  is  possible  to  predict  a  throat  coefficient  based 
on  the  correlation  (from  reference  (q)  ) 

Nu  pr  “■  =  ,021  Rep'® 

which  will  be  higher  than  predicted  on  the  flat  plate  correla¬ 
tion  (from  reference  (n)  ) 

-  1/3  8 

Nu  Pr  =  ,0295  Re_' 


For  this  case,  experimental  evidence  indicates  the  pipe  flow 
correlation  would  be  the  preferred  basis  for  predicting  the 
coefficients  in  the  throat  section.  This  is  not  to  be  con¬ 
strued  as  indicating  that  the  velocity  profile  in  the  throat 
section  of  the  nozzle  would  be  the  same  as  for  established 
pipe  flow.  Reference  (p)  reports  correlation  of  data  for 
rocket  combustion  chambers  and  nozzles  on  this  basis ,  It 
may  be  seen  in  Figs,  5a,  b  that  the  local  coefficients  ob¬ 
tained  for  the  throat  section  in  the  experiment  are  of  magni¬ 
tudes  corresponding  to  the  established  pipe  flow  behavior* 

The  method  of  testing,  using  small  thermal  plugs  to  cbtain 
the  local  heat  transfer  coefficients,  results  in  a  delayed  thermal 
boundary  layer  with  respect  to  the  origination  of  the  l:^drodynamic 
bouriary  layer.  This  condition  would  not  exist  in  a  nozzle  with 
heat  transfei  throughout.  Reference  (f)  gives  a  theoretically 
predicted  correction  for  this  effect  in  the  case  of  turbulent 
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boundary  layer  flow  over  a  flat  plate.  When  this  correction  is 
introduced,  the  experimental  results  agree  surprisingly  well  with 
the  flat  plate  equation 


Nu  Pr  =  .0295  Re  * 

X 

The  difference  between  Figs,  4a  and  4b,  5a  and  5b  is  the  intro¬ 
duction  of  this  compensating  factor. 

The  good  agreement  with  flat  plate  data  may  be  in  part 
fortuitous  since  the  correction  described  above  is  predicted 
for  flat  plate  flow  with  no  longitudinal  pressure  gradient; 
whereas,  tha  nozzle  flow  is  subjected  to  severe  axial  pressure 
gradients , 

A  boundary  layer  analysis,  utilizing  99^  of  the  free  stream 
velocity  as  the  thickness  of  the ^boundary  layer;  the  boundary 
layer  momentum  equation  of  von  Karmen  and  the  pressure  gradient 
of  nozzle  III,  has  been  made,  reference  (e) ,  The  results  indicate 
that  the  influence  of  nozzle  III  pressure  gradient  on  the  flat  wall 
hydrodynamic  boundary  layer  will  reduce  the  boundary  layer  thick¬ 
ness  in  the  throat  section  by  10%  when  compared  to  the  flat  plate 
hydrodynamic  boundary  layer  with  no  pressure  gradient.  This  re¬ 
sults  in  an  estimate  that  the  heat  transfer  coefficient  at  the 
throat  should  be  approximately  30^  higher  than  predicted  by  the  flat 
plate  correlation. 

From  l;hese  arguments,  it  is  possible  that  the  use  of  the 
thermal  plugs  tends  to  conceal  the  effect  of  the  pressure  gradient 
by  departing  from  the  correction  as  evaluated  from  reference  (f ) , 
This  remains  to  be  demonstrated. 

Until  more  experimental  evidence  is  obtained  it  is  re¬ 
commended  that  the  flat  plate  correlation 

Nu  Pr  "  4/3  ^  ^0295  Re^*® 

A 


be  used  to  prod  let  tho  luo,  I  heat  transfer  coefficients  in  the 
convergent  porMun  --i  a  do  Laval  nozzle. 

For  the  divtTKont  portion  of  the  de  Laval  nozzle,  re¬ 
ference  (d)  reports  correlation  of  results  for  a  particular 
nozzle  by  the  equation 

-  .2 

St  -  .0285  Re 
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The  length  in  the  Reynolds  Number  is  measured  from  the  throat 
of  the  nozzle • 

It  is  ^worthwhile  to  compare  the  results  of  reference  (d) 
with  the  results  of  this  experiment  since  the  two  sets  of  data 
overlap  at  the  throat  section  of  the  nozzles .  This  has  be^n 
done  using  the  local  hydraulic  diameter  as  the  characteristic 
length.  The  data  are  in  fair  agreement  as  can  be  seen  in  Fig.  6, 
Some  correction  of  data  of  reference  (d)  is  in  order  for  there 
appears  to  be  a  delayed  thermal  boundary  layer  in  this  test 
apparatus  also.  This  correction  would  bring  the  two  sets  of 
data  into  better  agreement. 
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EXPERIMENTAL  METHODS 


Modification  of  Experimental  Apparatus 

Results  reported  previously  were  for  cop^>er  plugs  buried 
in  the  curved  wall  only.  For  the  tests  reported  here  of  nozzles 
III,  IV,  and  V,  plugs  have  been  buried  along  the  centerline  of 
the  flat  wall  of  the  nozzles  as  well. 


Modification  of  Method  of  Reducing  Data 

The  procedure  used  in  reducing  the  data  reported  in  HS-2, 
reference  (b) ,  was  inadequate.  The  source  of  error  resulted  from 
the  heat  loss  from  the  copper  plug  to  the  bakelite  surrounding  it 
being  approximated  by  the  equation: 

12  "  V2 


where  h2  was  evaluated  by  taking  runs  with  no  flow  through  the 
nozzle.*^  This  value  was  then  subtracted  from  the  coefficient  ob¬ 
tained  by  reducing  the  transient  data,  for  flow  in  the  nozzle,  in 
the  standard  fashion. 


Actually  the  heat  loss  to  the  bakelite  is  represented  by 
the  equation: 

where  k^,  are  conditions  of  the  bakelite  at  the  copper  plug 

bakelite  interface,  and  for  the  actual  conditions  is  a 

variable  throughout  the  test  run.  ^  ° 


It  was  determined  by  a  graphical  solution  that  for  the  con¬ 
ditions  of  the  tests  and  with  the  materials  used,  \4.  reverses 

^dr^b 

the  time- 


sign  during  the  cooling  run  in  the  usable  portion  of 


temperature  cooling  curve .  At  the  time 


there  is  no 


heat  loss  to  the  bakelite  and  the  slope  of  the  time -log  temperatur 
curve  in  this  region  will  be 
evaluated , 


Vi 


from  which  h^  can  be 


V^( 


Use  was  made  of  this  fact  to  reduce  the  data  for  tnis  report. 
The  heat  transfer  coefficient  was  obtained  by  drawing  a  straight 
line  through  the  points  on  the  time-temperature  cooling  curve  in 


the  region  where  ,<^tN 

^Jr  • 

predicted  location  of 


was  predicted  to  be  near  zero.  The 


=  0 


was  justified  initially  by 


observing  that  the  time-temperature  curve  became  straight  in  this 
region,  which  it  should  for  no  heat  leakage.  An  example  time- 
temperature  curve  from  which  the  slope  has  been  defined  is  shown 
in  Fig .  (8) . 


Because  heating  time  is  important  in  predicting  the  point, 

0,  and  since  it  was  not  held  constant  in  previous  test 
b 

work,  the  data  available  was  not  recalculated.  Instead,  nozzle  III 
has  been  retested.  Nozzle  II  was  not  retested  since  it  was  decided 
that  nozzles  III,  IV,  and  V  would  give  adequately  descriptive  re¬ 
sults.  The  results  of  testing  nozzles  III,  IV,  and  V  are  reported 
here . 
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EXPERIMENTAL  RESULTS 


The  data  of  principal  interest  are  the  local  heat  transfer 
coefficients  at  the  throat  of  the  nozzles  and  the  influence  of 
the  various  wall  curvatures  on  these  coefficients.  Fig.  3  t)est 
describes  these  results.  From  this  figixpoe^it  may  be  observed 
that  for  the  high  mass  flows  (G  =  5. >2  x  10*"'  at  the  throat),  the 
local  coefficients  for  all  nozzles  for  both  the  contour  and  flat 
walls  just  after  the  throat  section,  may  be  described  by  h^  =  350  +10^. 
Note  however ,  that  the  ctifved  wall  coefi  icientf,  vary  over 
a  considerable  range,  h,  =  310  to  38C^  The  flat  wall  coefficients 
vary  only  from  ^  ^85  which  is  within  the  range  of  the 

■la 

experimental  accuracy  of  the  test  apparatus » 


The  flat  wall  of  the  nozzles  is  relatively  free  of  wall  curva 
tiire  6ff«»cts,  and  hence  these  results  Indicate  that  the  variation 
of  negatj:.ve  pressure  gtadients  of  the  magnitude  found  in  the  test 
nozzles  has  only  minor  influence  on  the  local  heat  transfer  coeffi 
cients  at  the  throat. 


The  effect  of  wall  curvature  is  more  significant.  Since  there 
is  no  difference  in  wall  curvature  in  the  three  nozzles  precisely 
at  the  throat  or  thereafter,  the  variation  in  the  local  coefficient 
(310  to  380)  among  the  three  nozzles  is  assumed  to  be  the  relative 
effect  of  different  curvatures  upstream  in  the  nozzles  that  has  per¬ 
sisted  to  the  throat. 


Greater  evidence  of  the  influence  of  wall  curvature  can  be  ob¬ 
served  by  examining  the  curved  wall  coefficients  of  the  nozzles  up¬ 
stream  of  the  nozzle  throats.  Here  it  is  seen  that  the  wall  curva¬ 
ture  has  a  marked  effect  on  the  local  coefficient,  A  concave  curva¬ 
ture  such  as  in  nozzle  IV  produces  a  coefficient  k0%  lower  than  that 
realized  by  the  flat  wall  in  the  same  location  along  the  nozzle  axis. 
The  converse  is  true  for  the  sharp  convex  curvature  found  in  nozzle  V, 
where  the  curved  wall  coefficient  is  about  k0%  higher  than  that  real¬ 
ized  by  the  flat  wall  in  the  same  location, 

FigSi,  ^a,  b,  and  5a,  b  demonstrate  that  the  heat  transfer 
behavior  In  the  convergent  portion  of  the  nozzles  more  nearly  ap¬ 
proximates  flat  plate  rather  than  pipe  flow  behavior.  The  results 
in  Figs,  ba  and  5a  are  high.  This  is  due  to  the  fact  that  for  the 
experimental  tests  small  thermal  plugs  were  used  which  results  in 
a  delayed  thermal  boundary  layer  with  respect  to  the  initiation  of 
the  hydrodynamic  boundary  layer.  This  effect  has  been  taken  into 
account  in  plotting  the  data  in  F.igs,  ^-b  and  5b.  The  results  in 
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Fig.  ^-b  show  good  agreement  with  the  flat  plate  correlation  for 
heat  transfer 

-  1/3  8 

Nu  Pr  =  .0295  Re^*'^ 

Fig.  6b  demonstrates  that  prediction  of  heat  transfer  coefficients 
in  the  convergent  portion  of  a  nozzle  using  the  conventional  estab¬ 
lished  pipe  flow  correlation  would  lead  to  results  considerably  below 
the  actual  conditions.  However,  it  is  to  be  noted  that  at  the  throat 
section  the  coefficients  agree  with  established  pipe-flow  coefficients. 

The  experimental  results  obtained  for  the  throat  section  of  the 
nozzles  have  been  compared  to  those  reported  in  reference  (d)  for 
the  throat  and  divergent  portion  of  a  particular  nozzle .  This  com¬ 
parison  is  shown  in  Fig .  6  using  the  local  hydraulic  diameter  as 
the  characteristic  length  in  evaluating  the  Reynolds  Number .  The 
experimental  data  of  this  report  have  been  modified  by  the  correc¬ 
tion  factor  for  the  delayed  thermal  boundary  layer,  as  described 
earlier.  The  data  of  reference  (d)  have  been  altered  by  the  author, 
in  that  the  viscosity  has  been  evaluated  on  an  estimated  average 
value  between  the  wall  temperature  and  the’  stream  temperature. 

The  original  data  was  reported  using  the  viscosity  based  on  the 
stream  stagnation  temperature.  These  data  are  in  fair  agreement. 

The  test  apparatus  used  in  reference  (d)  appears  to  introduce 
a  delayed  thermal  boundary  layer.  Compensation  for  this  effect 
would  bring  the  data  into  better  agreement . 

In  Fig.  7  the  heat  transfer  coefficients  for  the  plugs  lo¬ 
cated  on  either  side  of  the  throat  section  are  plotted  using 
the  conventional  flat  plate  dimensionless  groups  and  introducing 
the  correction  factor  just  described. 

Fig.  10a  compares,  for  a  particular  mass  flow,  the  experi¬ 
mental  local  heat  transfer  coefficients  obtained  along  the  flat 
wall  of  nozzle  III  to: 

a)  The  local  heat  transfer  coefficient  predicted  by  the 
flat  plate  correlation 


Nu^r"  =  .0295  Re^*® 

b)  The  local  heat  transfer  coefficient  as  predicted  from 
reference  (ra)  for  unestablished  flow  in  pipes , 
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c)  The  local  heat  transfer  coefficient  predicted  from  the 
boundary  layer  thickness  as  influenced  by  the  pressure 
gradient  of  nozzle  III. 

d)  The  experimental  local  heat  transfer  coefficient  cor¬ 
rected  as  suggested  in  reference  (f)  for  the  delayed 
thermal  boundary  layer  effect. 

Figs .  10b  and  c  compare  the  local  predicted  flat  plate  coeffi¬ 
cients.,  the  experimental  flat  wall  coefficients,  and  the  experi¬ 
mental  local  coefficients  taking  into  account  the  delayed  thermal 
boundary  layer  effect.  These  comparisons  are  for  the  flat  "-falls 
of  nozzles  IV^and  V  respectively,  and  the  same  mass  flows 
(G  =  5.2  -x  10°  at  the  throat).  It  is  seen  that  in  all  cases,  in 
the  throat  region,  the  corrected  predicted  results  agree  sur¬ 
prisingly  well  with  the  experimental  results . 
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Figure  la 


Pressure  ratio  versus  distance  from  nozzle 
entrance  for  test  nozzle  contours  III,  iv,  v. 


Figure  lb 


Nozzle  area  versus  distance  from  nozzle  entrance 
for  test  nozzle  contours  III,  IV,  V. 


Figure  2a 


Picture  of  bakelite  contour  with  embedded 
copper  plugs  and  brass  counterpart,  nozzle  IV, 


Figure  2b 


Picture  of  nozzle  IV  with  one  flat  wall  removed. 
Note  copper  plugs  buried  in  flat  wall. 
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PRESSURE  RATIO  if/Po)  CROSS  SECTIONAL  AREA  OF  NOZZLE  (IN.*)  •  1/2  NOZZLE  WIDTH 


DISTANCE  FROM  NOZZLE  ENTRANCE  (INCHES)  FIGURE  I  b 


n.G. 


riGZb 


Figure  3  Local  apparent  heat  transfer  coefficients  versus 
distance  from  nozzle  entrance  for  test  nozzles 
III,  IV,  V. 

O  Points  are  coefficients  obtained  from  contoured 
walls  of  nozzles. 

tZl  Points  are  coefficients  obtained  from  flat 
walls  of  nozzles, 

-  -  Curves  are  for  high  mass  flow. 

G  =  5.2  X  10^  at  nozzle  throats  and  M  =  1  , 

-  Curves  are  for  low  mass  flow, 

G  =  3,5  X  10''^  at  nozzle  throats  for  nozzles  IV 
and  V  and  M  =  1. 

G  =  2.6  X  10^  and  M  =  ,85  at  nozzle  throat  for 
nozzle  III, 


LOCAL  HEAT  TRANSFER  COEF.  BTU/HR. 


DISTANCE  FROM  NOZZLE  ENTRANCE  INCHES  FIGURE  3 


Figure 


Figure 


Figure 


Figure 


4a  Experimental  data  for  nozzles  III,  IV,  V  are 

plotted  in  dimensionless  form,  using  distance 
measured  from  the  nozzle  entrance  along  nozzle 
profile  as  characteristic  length  dimension. 

Nu  Pr""  Versus  Re 


4b  Same  as  Figure  4a  except  that  compensating 
factor  for  delayed  thermal  boundary  layer 
(from  reference  f)  has  been  used  to  correct 
data  for  this  effect. 


5a  Experimental  data  for  nozzles  III,  IV,  V  are 
plotted  in  dimensionless  form  using  local  hy¬ 
draulic  diameter  as  the  characteristic  length 
dimension. 

Nu  Pr”  Versus  Re^ 


5b  Same  as  Figure  5a  except  that  the  compensating 
factor  for  delayed  thermal  boundary  layer 
(from  reference  f)  has  been  used  to  correct 
data  for  this  effect. 
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Figure  b 


Figure  7 


Figure  8 


Data  for  the  throat  sections  of  test  nozzles  III, 
IV.  V,  are  plotted  in  dimensionless  form  after 
using  the  compensating  factor  for  the  delayed 
thermal  boundary  layer,  reference  (f).  Also  the 
data  from  reference  (d;  for  the  throat  section 
and  divergent  section  of  a  nozzle  are  plotted  on 
the  same  dimensionless  basis.  Sections  1,  2^  3, 

4,  of  reference  (d)  refer  to  the  throat  section  and 
the  subsequent  three  sections  comprising  the  di¬ 
vergent  portion  of  the  nozzle  respectively. 

Data  for  two  plugs  nearest  throat  section  on  both 
contour  and  flat  walls  of  nozzles  III,  IV.  V  are 
plotted  on  dimensionless  basis  using  profile  length 
from  nozzle  entrance  as  characteristic  length. 

The  data  have  been  compensated  for  the  delayed 
thermal  boundary  layer  effect. 

Sample  Temperature -Time  cooling  curve  of  copper 
plug  from  which  the  heat  transfer  coefficient  is 
deduced. 

*■  versus  cooling  time, 

p  &w 
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Figure  9.  Curves  are  for  a  particular  mass  flow  through 

nozzle  III  and  are  plotted  versus  distance  from 
nozzle  entrance, 

(A)  Flat  plate  turbulent  boundary  layer  thickness. 

(B)  Flat  plate  boundary  layer  thickness  as  in¬ 
fluenced  by  pressure  gradient  of  nozzle 

reference  (e)  . 

(C)  Mass  flow  velocity. 

Figure  10a, 

Curves y  for  a  particular  mass  flow  through 

nozzle  III ,  pertain  to  flat  wall  of  the  nozzle 
and  are  plotted  versus  distance  from  nozzle 
entrance . 

(D)  Experimental  local  heat  transfer  coefficients. 

(E)  Predicted  local  coefficients  on  flat  plate 
basis  after  introducing  the  pressure  gradient 
of  nozzle  III,  reference  (e) , 

(F)  Predicted  local  coefficients  from  flat  plate 
correlation. 

(G)  Local  coefficients  predicted  for  unestablished 
flow  in  pipes ,  reference  (m) , 

•  Points  are  experimental  coefficients  cor¬ 
rected  for  delayed  thermal  boundary  layer 
effect,  reference  (f)  , 


Figures  10b,  c. 

Flat  wall  heat  transfer  coefficients  are  plotted 
versus  distance  from  nozzle  entrance  for  nozzles 
IV,  V,  respectively,  and  for  c 

G  =  5.2  X  10^ 


at  nozzle  ttooats , 


O  Points  are  experimental  data. 


•  Points  are  experimental  data  corrected  for 
influence  of  delayed  thermal  boundary  layer. 
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DISTANCE  FROM  NOZZLE  ENTRANCE  FEET  FIGURE  10  b  DISTANCE  FROM  NOZZLE  ENTRANCE  ALONG  FUT  VSlL  FEET  FIQURE  IQc 


DISCUSSION 


Before  discussing  the  results  as  they  would  apply  to  a  con¬ 
ventional  nozzle  where  heat  transfer  takes  place  continuously 
along  the  nozzle  wall,  it  is  necessary  to  give  some  consideration 
to  the  conditio'  that  in  these  tests  the  hydrodynamic  and  thermal 
boundary  layers  did  not  originate  simultaneously.  The  hydrodynamic 
boundary  layer  starts  at  the  entrance  to  the  nozzle  whereas  the 
thermal  boundary  layer  starts  at  the  upstream  edge  of  the  heated 
plug.  Neglecting  for  the  moment  the  effects  of  pressure  gradient, 
mass  acceleration,  and  wall  curvature j  it  is  apparent  that  due  to 
the  delayed  thermal  boundary  layer  the  heat  transfer  coefficients 
obtained  by  this  experimental  apparatus  will  be  than  would 

be  expected  in  the  duplicate  nozzle  with  heat  tj  '  :  '^ughout , 

Then  the  problem  arises,  to  estimate  quantitatively  *ie  eiiect  of 
the  delayed  thermal  boundary  layer  on  the  local  heat  transfer 
coefficient , 

However,  before  this  problem  can  be  resolved,  it  is  necessary 
to  have  some  basis  for  predicting  the  local  heat  transfer  coeffi¬ 
cient  for  heat  transfer  throughout  the  nozzle.  Two  bases  of  pre¬ 
diction  are  possible.  One,  that  the  nozzle  is  so  short,  L/D  of 
approximately  10  based  on  the  throat  hydraulic  diameter,  that  the 
flow  is  unestablished  and  hence  is  essentially  the  same  as  flow 
over  a  flat  plate.  The  other,  that  the  flow  is  nearly  established 
and  hence  the  pipe  flow  correlation  is  a  better  basis  of  analysis , 

That  the  flat  plate  basis  is  a  superior  starting  point  is 
apparent  by  examining  Figs,  b,  and  5a,  b.  The  data  shows 
marked  evidence  toward  correlation  using  distance  from  the  nozzle 
entrance  as  the  characteristic  length  whereas  there  is  no  correla¬ 
tion  when  using  the  local  hydraulic  diameter  as  the  characteristic 
length. 

Then  the  flat  plate  correlation  suggested  by  Colburn, 
reference  (n) , 


St  Pr^^^  -  .0295  Re„"*^ 

Jv 


will  be  utilized  in  the  effort  to  estimate  the  influence  of 
the  delayed  thermal  boundary  layer  on  the  local  heat  transfer 
coefficient , 
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The  basis  of  analysis  could  also  be  th'^  pipe  entrance 
flow  behavior  as  predicted  by  reference  (m)o  However,  this 
equation  is  considerably  more  cumbersome  to  work  with  and  gives 
results  that  agree  quite  well  with  the  flat  plate  predictions 
(see  Figo  10a )o  Therefore,  the  flat  plate  correlation  will  be 
usedo 


Examination  of  Figo  4  reveals  that  the  experimental  results 
are  higher  by  100^  for  the  plugs  near  the  nozzle  entrance  and 
about  60Jf  higher  for  the  plugs  near  the  throat  section  when  com¬ 
pared  to  the  flat  plate  correlationo 

The  plugs  near  the  nozzle  entrance  are  not  influenced  by 
pressure  gradient  or  wall  curvature o  Also,  they  are  near  the 
I’nitiation  of  the  hydrodynamic  boundary  layer  so  the  flat  plate 
flow  behavior  should  be  valid »  Then  the  remaining  faptors  caus¬ 
ing  these  data  to  be  100^  high  are  the  delayed  thermal  boundary 
layer  and  abnormal  nozzle  entrance  conditions o  The  nozzle  en¬ 
trance  from  the  static  chamber  was  made  with  large  radius  contours <> 
Furthermore,  these  contours  were  varied  and  produced  no  signi¬ 
ficant  change  in  the  experimental  results »  Therefore,  it  is 
indicated  that  the  delayed  thermal  boundary  layer  is  x;ha  princ.  ’.pal 
cause  for  the  high  data*, 

Consider  the  two  possible  extreme  combinations  of  thermal 
and  hydrodynamic  boundary  layer  origination  with  respect  to  the 
thermal  plugo  The  extreme  which  would  give  the  lowest  heat  trans¬ 
fer  coefficient  for  the  plug  would  be  if  both  thermal  and  hydro- 
dynamic  boundary  layers  originated  at  the  nozzle  entrance „  The 
alternate  extreme  or  one  yielding  the  highest  local  heat  transfer 
coefficient  would  be  if  both  the  thermal  and  hydrodynamic  boundary 
layers  originated  at  the  leading  edge  of  the  plugo 

Based  on  flat  plate  behavior,  for  all  stream  properties 
constant  and  free  stream  velocities  identical,  the  local  heat 
transfer  coefficient  h^^  is  proportional  to  o2  for 

turbulent  boundary  layer o 

For  the  plugs  near  the  nozzle  entrance  this  gives  a  varia¬ 
tion  of  h^,  based  on  the  two  extremes  just  described,  of  approx¬ 
imately  lo??  the  latter  coefficients  being  70^highero  Then. 
if  the  actual  plug  behavior  were  near  the  upper  limit  of  the  wo 
extremes  it  is  expected  that  the  data  would  be  approximately 
70^  higher  than  the  Colburn  correlation  predicts o  That  the  re¬ 
sults  for  the  first  plugs  are  on  an  average  30^  higher  than  this 
upper  limit  cannot  at  this  time  be  accounted  foro 
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That  the  coefficients  obtained  experimentally  should  be 
nearer  the  upper  limit  than  the  lower  limit  can  be  justified* 

The  more  significant  of  the  two  boundary  layer  thicknesses  with 
respect  to  heat  tran*^fer  is  the  thermal  boundary  layer*  Of  the 
two  extremes  considered,  the  one  with  the  thermal  boundary  layer 
starting  at  the  leading  edge  of  the  thermal  plug  is  much  nearer 
the  actual  test  conditions* 

This  argument  was  further  substantiated  by  placin,^  thermal 
plug,  identical  with  the  test  plugs,  in  the  wall  of  a  c..rcular 
pipe*  Here  the  flow  was  establishea  turbulent  flow*  ^he  two 
extremes  between  which  the  plug  coefficient  should  lie  in  this 
case  were§  for  the  lower  limit,  the  heat  transfer  coefficient 
predicted  from  the  turbulent  pipe  flow  correlation  and  for  the 
upper  limit,  the  coefficient  assuming  both  the  thermal  and  hy¬ 
drodynamic  boundary  layers  started  at  the  leading  edge  of  the  plug 
From  the  previous  argument  the  results  should  be  mar  the  upper 
extreme  value*  This  is  the  case*  The  upper  value  was  predicted 
as  250< higher  than  the  lower  limit  and  the  experimental  values 
were  230^  higher  than  this  lower  limit*  These  results  agree  with 
theory*  Reference  (m)  predicts  that  the  data  should  be  230JI? 
above  the  lower  limit  due  to  the  effect  of  unestablished  thermal 
boundary  layer* 

From  this  analysis  it  is  concluded  that  the  maximum  delayed 
thermal  boundary  layer  effect  on  the  experimental  values  for  the 
plugs  near  the  nozzles’  entrance  would  be  to  raise  them  approxi¬ 
mately  60“70  ^higher  than  predicted  by  flat  plate  data* 

In  reference  (f)  a  theoretical  prediction  is  made  of  the  ef¬ 
fect  of  the  delayed  thermal  boundary  layer  (wall  temperature  dis¬ 
continuity)  for  turbulent  flow  over  a  flat  plate*  The  results 
predicted  for  the  plugs  near  the  nozzle  entrance,  based  on  this 
reference,  should  be  40-60 higher  than  those  for  the  case  of 
heat  transfer  throughout  the  nozzle*  For  the  plugs  near  the 
throat  of  the  nozzle  this  reference  predicts  that  the  results 
will  be  60-80^ higher  than  for  the  case  with  heat  transfer  through 
out*  These  values  are  in  agreement  with  the  foregoing  analysis 
and  have  been  used  to  eliminate  the  effect  of  the  delayed  thermal 
boundary  layer  in  comparing  the  data  to  correlations  where  the 
thermal  and  hydrodynamic  boundary  layers  originate  simultaneously* 

The  plugs  farther  removed  from  the  nozzle  entrance  are  under 
the  influence  of  pressure  gradient,  mass  acceleration  and  wall 
curvature*  That  they  show  a  strong  tendency  to  c orrelate  is  in¬ 
dicative  that  these  variables  do  not  exert  extreme  infa.aences 
with  respect  to  the  local  heat  transfer  coefficients* 
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The  delayed  thermal  boundary  layer,  based  on  the  foregoing 
discuLsion,  should  cause  these  results  also  to  be  high  by  an 
approximate  factor  of  60“80^o  This  correction  agrees  with  the 
experimental  values o  Hov/ever,  it  is  known  that  the  negative 
pressure  gradient  has  some  influence  on  reducing  the  hydro- 
dynamic  boundary  layer  thickness  and  hence  should  raise  the  local 
heat  transfer  coefficiento 

The  effect  of  the  negative  pressure  gradient  has  been  quan¬ 
titatively  approximated  for  nozzle  III  for  a  particular  mass  flow. 
The  von  Karman  boundary  layer  momentum  equation,  the  l/7th  power 
velocity  distribution, and  the  pressure  gradient  of  nozzle  III  were 
utilized  and  integrated  graphically  to  predict  the  boundary  layer 
thickness  as  a  function  of  distance  from  the  nozzle  entrance,  ref¬ 
erence  (e)«  The  results  are  shown  in  Fig*  9  compared  to  the 
boundary  layer  thickness  based  on  flat  plate  flow,  turbulent  bound¬ 
ary  layer  and  no  pressure  gradiento  It  is  to  be  seen  that  the 
pressure  gradient  influence  reduces  the  boundarv  layer  thickness 
at  the  throat  of  the  nozzle  to  approximately  l/+  of  the  value  for 
no  pressure  gradientc 

Using  the  boundary  layer  thickness  predicted  by  this  analysis, 
the  local  heat  transfer  coefficient  was  evaluated  by  the  relation 
derived  by  Prandtl,  which  predicts  h,  as  a  function  of  fluid  prop¬ 
erties,  wall  shear  stress,  maximum  velocity  in  the 

laminar  sublayer,  and  free  stream  velocityo 

The  results  are  comparei  in  Figo  10  to  the  results  predicted 
by  the  flat  plate  correlation  and  the  results  obtained  by  experi¬ 
ment  for  the  flat  side  of  nozzle  IIIo  It  is  seen  that  the  pre¬ 
dicted  coefficient,  considering  pressure  gradient,  is  approximately 
higher  in  the  throat  region  than  the  coefficient  based  on 
flat  plate  basis  slone 

Now  if  the  delayed  thermal  boundary  factor  remained 

as  predicted  before  consideration  of  pressure  gradient  and  the 
pressure  gradient  influenced  the  coefficient  as  just  described, 
the  experimental  coefficients  for  the  plugs  near  the  throat  should 
be  approximately  100 ^higher  than  predicted  on  the  flat  plate 
basis  onlyo  However,  tho  delayed  thermal  boundary  layer  effect, 
since  it  is  a  function  of  the  thermal  and  hydrodynamic  boundary 
layers,  will  also  be  influenced  by  the  negative  pressure  gradient* 
It  should  be  reduced  as  a  result  of  the  reduction  in  thickness 
of  the  hydrodynamic  boundary  layer  thickness*  An  estimate  of 
this  change j  based  on  the  graphically  predicted  boundary  layer 
thickness,  is  that  the  experimental  results  should  be  20-30^ 
higher  than  for  a  continuously  heated  plate*  Combining  the  two 
influences,  the  delayed  thermal  boundary  layer  and  the  pressiire 
gradient,  indicates  that  the  results  for  the  plugs  near  the 
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throat  sectio'n  should  be  50~60%  higher  than  predicted  by  the 
flat  plate  correlationo  This  is  in  agreement  with  the  results 
obtained  experimentallyo 

It  is  to  be  observed  from  the  preceding  discussion  that 
it  is  possible  for  the  delayed  thermal  boundary  layer  factor 
to  conceal  the  effect  of  the  negative  pressure  gradient  by  de¬ 
creasing  from  the  prf  dieted  value  in  amount  equal  to  the  increase 
in  the  local  coefficient  due  to  the  pressure  gradient  effect o 
Whether  or  not  this  is  so  remains  to  be  proveno 

At  this  time,  available  theory  does  not  permit  the  pre¬ 
diction  of  the  local  heat  transfer  coefficient  under  the  in¬ 
fluence  of  both  pressure  gradient  and  wall  curvature »  From  the 
experimental  results,  see  Figo  3,  it  may  be  seen  that  in  no 
case  did  the  wall  curvature  produce  a  departure  between  the 
flat  and  contoured  walls  of  more  than  40^o  Also  when  the  re¬ 
sults  from  the  two  walls  are  plotted  in  dimensionless  form  as 
in  FigSo  4a,  b  and  5a,  b,  it  is  to  be  seen  that  they  plot  in 
the  same  manner  although  the  scatter  of  the  results  for  the  con¬ 
toured  walls  is  greatero 

It  was  anticipated  initially  that  the  results  for  the 
nozzles,  particularly  near  the  entrance,  would  be  laminar, 
since  the  length  Reynolds  NOc  was  less  than 

5  X  10^0  In  no  case  were  laminar  heat  transfer  coefficients 
obtainedo  This  is  attributable  to  three  possible  causes o 
One,  that  the  boundary  layer  is  being  heated  and  hence  tends 
toward  early  transition^  two,  that  even  though  the  plugs  are 
machined  flush  with  the  bakelite,  the  joint  between  the  bake- 
lite  and  the  copper  plug  is  not  hydrodynamically  clean  and 
hence  tends  t^  promote  turbulence  in  the  boundary  layer;  three, 
that  the  stream  entering  the  nozzle  has  sufficient  degree  of 
turbulence  to  promote  early  transitiono 

The  topic  of  established  versus  non-established  flow  in 
the  convergent  portion  of  the  de  Laval  nozzles  should  be  dis- 
cussedo  Taking  into  account  the  correction  factors  discussed 
previously  it  is  to  be  observed  in  Flgc  5  that  the  data  in  the 
throat  region  are  in  agreement  with  the  coefficients  predicted 
for  established  flow  in  pipes*  It  has  teen  experimentally 
determined,  references  (1)  and  (m).  that  at  an  L/D  of  approxi¬ 
mately  10  for  the  test  flow  conditions,  the  heat  transfer  coeffi¬ 
cient  is  within  lO^iof  that  predicted  for  established  pipe  flow* 

If  the  nozzle  throat  hydraulic  diameter  is  used,  the  maximum 
L/D  is  approximately  lO*  Therefore,  the  indication  that  tlie 
coefficients  c.t  the  throat  are  approaching  established  pipe 
flow  coefficients  is  to  be  expected*  Therefore,  if  throat 
coefficients  are  of  interest  to  a  nozzle  designer,  the  L/D 
value  of  the  throat  should  be  a  design  parameter* 
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Reference  (d)  has  reported  results  for  the  throat  and  di¬ 
vergent  sections  of  a  particular  de  Laval  nozzle.  Based  on  the 
author’s  interpretation  of  an  appropriate  Reynolds  number  (using 


t 


stream 


^wall 


2 


to  evaluate  viscosity) 


these  results  are  shown  in  Fig,  6,  Section  1  data  of  reference 
(d)  are  average  coefficients  in  the  nozzle  throat  section  and 
hence  can  be  compared  with  the  throat  section  data  reported  here. 
It  is  to  be  seen  that  these  data  are  in  fair  agreement. 

In  summary,  it  is  anticipated  that  the  local  heat 
transfer  coefficients  in  the  convergent  portion  of  the  nozzle 
will  agree  essentially  with  flat  plate  coefficients,  Hiowever, 
due  to  curvature  and  pressure  gradiej  ^  effects  they  may  vary 

-  30^^  from  this  value  as  the  nozzle  \  oat  is  approached.  In 
the  throat  region  they  may  be  igh  as  a  result  of  the 

negative  pressure  gradient  although  this  is  not  proven.  Further¬ 
more,  it  is  recommended  that  if  calculation  on  a  pipe  Reynolds 
Number  basis  gives  a  higher  coefficient  than  the  flat  plate 
Reynolds  Number  basis,  it  should  be  used  for  predicting  the  throat 
heat  transfer  coefficient,  as  a  "conservative "  design  procedure. 

The  proposed  future  wor :  on  this  project  becomes  apparent 
from  the  discussion  of  results  obtained  thus  far* 

a)  Discover  why  the  nozzle  entrance  coefficients  are 
higher  by  30/l'  than  can  be  justified, 

b)  Isolate  the  influence  of  pressure  gradient  on  the  delayed 
thermal  boundary  layer  so  that  the  data  reported  can  be 
unquestionably  used  for  the  conditions  of  continuous 
heat  trn^sfer  in  the  nozzle. 
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APPENDIX 


METHOD  OF  DETERMINING  THE  HEAT  TRANSFER  COEFFICIENT 
FROM  THE  TIME -TEMPERATURE  HISTORY  OF  THE  THERMAL  PLUG 


The  method  of  predicting  the  heat  transfer  coefficient  to 
the  gas  stream  using  the  thermal  plug  and  transient  technique 
would  be  a  straight  forward  matter,  if  during  the  heating  and 
the  cooling  of  the  plug,  the  heat  loss  to  the  bakelite  matrix 
holding  the  plug  were  negligible. 


For  this  transient  case  during  cooling  the 
coefficient  can  be  evaluated  from  the  equation. 


loggC 


^9 

^9"  -  '*^aw 


0”  -  0* 


heat  transfer 


(1) 


By  plotting  (t^-  t  )  versus  0  on  semilog  paper  and  using  the 
slope  of  this  ^  ^  plot  the  heat  transfer  coefficient  h-, 

can  be  evaluated  directly. 

This  is  not  the  case  in  the  experimental  apparatus  under 
the  conditions  of  test.  During  the  cooling  of  the  plugs  the 
energy  absorbed  by  the  bakelite  amounts  to  approximately  20% 
of  the  total  thermal  energy  stored  in  the  plugs  during  the  heating 
process . 

The  value  of  20%  was  predicted  by  a  graphical  solution  of 
the  time -temperature  history  in  the  bakelite.  The  graphical 
technique  used,  is  presented  in  reference  (a)  and  is  a  modifica¬ 
tion  for  radial  heat  flov;  of  the  Schmidt  method. 

Furthermore,  if  the  heat  loss  from  the  plug  to  the  phenolic 
were  a  constant  with  respect  to  time  or  varied  with  respect  to 
(tg-  t  )  some  modification  of  equation  (1)  might  be  usable.  But 
again,  rhis  is  not  the  case,  Therefor3  some  other  technique  is 
necessary  to  evaluate  the  heat  transfer  coefficient  to  the  gas 
stream. 

In  the  experimental  apparatus,  negleccing  heat  transfer 
from  the  back  side  of  the  plugs  which  was  found  to  be  in  all  cases 
less  than  3%  of  the  heat  transfer  rate  to  the  gas  stream,  heat 
flows  from  the  plug  via  two  paths .  One  into  the  high  velocity 
gas  stream,  represented  by  the  equation: 
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(2) 


The  other  is  into  the  bakelite  matrix  surrounding  the  plug,  rep¬ 
resented  by  the  equation: 


qa 


/<^t, 
^2^5r‘ 


t 

Where  is  the  radial  temperature  gradient  in  the  bakelite 

at  the  contact  face  between  the  bakelite  and  the  copper  plug. 

Using  the  graphical  solution  mentioned  previously^  and  the 
recorded  temperature  history  of  the  plug,  it  can  be  shown  that 
for  the  heating  time  involved  in  this  experiment  q^  -  qp  at  the 
start  of  the  plug  cooling  or  test  run.  This  further  amplifies 
the  need  for  taking  into  account  correctly  the  heat  loss  to  the 
bakelite  if  the  gas  heat  transfer  coefficient  is  to  be  evaluated 
with  reasonable  accuracy. 


Now  note  the  quantities  variable  with  time  in  equations  (2) 
and  (3).  In  (2)  (tg-  t  )  decreases  with  time  of  cooling j  likewise, 
in  (3)  (<)t/cir),  decreasss  with  time  of  cooling  for  the  test  condi¬ 
tions,  It  is  possible  to  write  (c}t/<)r).  as  t^-  t  where  (tr,-  t^) 
is  the  variable.  This  presupposes  °  ^  ^  ” 

that  the  plug,  bakelite  interface 
thermal  resistance  is  negligible  andAr<<r. 


Ar 


Then  it  remains  to  compare  (t^-  t„)  and  (t^  -  t  )  for  the 
cooling  time  of  the  plug.  Due  to  ^  ^  the  ^  ^  thermal 

energy  storage  in  the  bakelite,  t  will  increase  with  time  whereas 
^aw  remain  constant.  Therefore,  as  the  plug  cools  (tp  -  t  ) 

will  approach  zero  more  rapidly  than  (tp-  t  )  since  ^  ^ 

t  must  of  necessity  be  greater  than  t_^,  at^the^oeginning  of  the 
^  cooling  time. 


Also,  since  (tp-  t  )  is  not  zero  when  tp^  t  ,  thermal  energy 
will  continue  to  flow  from  the  plug  into  the  gas  stream  with  the 
result  that  (tp  -  t  )  becomes  negative.  Then  it  is  seen  that 
heat  flow  rate  from-^the  plug  during  the  initial  part  of  the  cooling 
curve  is  given  by  the  equation  q™  -  q.  ’»*  q^  j  whereas  during  the 
■"atter  portion  it  will  be  q^  =  -  qp. 


When  tp-  t  ,  qp  =0  and  equation  (1)  is  valid.  This  sugr 
that  if  the  time  0  at  which  qp  -  0  could  be  predicted  precis- 
and  were  in  a  suitable  portion  of  the  cooling  curve,  equatic 
could  be  used  at  this  point  on  the  cooling  curve  to  predict 
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However,  it  is  only  possible  to  predict  this  point  approximately 
since  a  graphical  solution  is  the  only  means  of  prediction.  But 
q_2  is  approximately  zero  compared  to  q-,  for  a  considerable  time 
interval  on  either  side  of  the  point  qi  =  0.  Specifically,  for 
qp  -  10^  q^  ^0  ^  10  sec.  for  this  experimental  apparatus.  Since  the 
total  usable  portion  of  the  cooling  run  is  30  to  60  seconds  dura¬ 
tion,  it  is  seen  that  A©  -10  sec,  covers  a  considerable  portion 
of  this  curve.  Also,  from  graphical  p>'edicticn  this  range  of 
interest  is  in  the  usable  portion  of  the  curve ,  This  is  shown 
in  Figo  8  which  is  a  typical  time -temperature  cooling  curve. 

Therefore,  by  confining  the  analysis  of  the  time -temperature 
cooling  curve  to  the  region  just  specified,  it  is  possible  to 
use  equation  (1)  and  evaluate  the  gas  stream  heat  transfer  coef¬ 
ficient  directly. 

Before  this  can  be  done  it  is  necessary  to  predict  in  what 
portion  of  the  curve  qp  ^  10^  q, ,  This  will  depend  on  the 
following  variables  s 

(a)  Plug  heating  time 

(b)  Maximum  temperature  of  the  plug  above 

(c)  Gas  film  heat  transfer  coefficient 

(d)  Ratio  of  area  of  plug  exposed  to  gas  stream  and  bakelite 

By  appropriate  experimental  technique  (a) ,  (b) ,  (d)  can  be 

held  constant  for  all  tests  and  plugs Therefore  the  location 

on  the  cooling  curve  of  ^2  ”  Qn  will  be  dependent  only  on 

hn .  From  a  series  of  graphical  predictions  for  various  values 

Oi  ^  the  range  of  interest,  the  location  on  the  cooling  curve 

of  qp  0  was  determined.  It  was  found  that  for  the  range  of 

h-,  or  concern  the  point  of  qp  ^  0  w<,s  approximately  a  fixed  uortion 

or  the  total  temperature  rise  of  the  plug  above  t  . 

aw 

This  basis  of  interpreting  the  data  was  then  used  to  evaluate 
sloug  with  physical  observation  of  the  time -temperature  cooling 
carve  to  see  which  group  of  points,  in  the  region  q^  -  0,  would 
yield  the  best  straight  line,  ^ 

Fig.  (11)  shows  the  graphical  solution  for  the  predicted  time- 
temperature  history  during  the  heating  of  a  plug . 

Fig,  (12)  is  the  graphical  solution  for  the  cooling  time 
tei^erature  history.  Note  the  point  of  qp  =■  0  is  at  approximately 
^^^tg  -  ol5, 

max 

As  a  further  check  on  the  validity  rf  the  graphical  solution, 
thermocouples  were  buried  in  the  bakelite  at  specified  radial 
distances  from  one  of  the  test  plugs.  Thus  time -temperature 
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histories  of  these  radial  points  were  also  obtained.  Superimposed 
on  these  results  were  the  results  predicted  by  the  graphical 
solution.  The  comparison  is  made  in  Fig.  (13)  • 

The  thermocouples  were  buried  near  the  gas  stream  face  of  the 
bakelite  to  determine  whether  or  not  the  heat  flow  from  the  plug 
into  the  bakelite  was  essentially  radial.  It  is  to  be  observed 
that  the  graphically  predicted  time -temperature  history  compares 
favorably  with  the  measured  results  for  the  thermocouple  located 
l/l6  inch  from  the  test  plug.  Whereas,  the  measured  results  are 
appreciably  lower  than  the  predicted  for  the  thermocoup] e  l/8 
inch  from  the  plug ,  This  is  no  doubt  due  to  the  heat  loss  from 
the  bakelite  to  the  gas  stream. 

From  these  results  it  may  be  concluded  that  the  heat  flow 
from  the  plug  to  the  bakelite  is  essentially  radial  near  the 
plug  and  that  the  graphical  solution  is  satisfactory  for  these 
test  conditions , 


Figure  11  Copper  plug  temperpture  rise  versus  radial  distance 
on  a  log  scale  w?th  time  as  a  parameter. 


Figure 


Figure 


Graphical  solution  of  time -temperature  history  of 
bakelite  surrounding  a  copper  plug  while  the  plug 
is  being  heated. 


12  Copper  plug  temperature  decline  versus  radial 

distance  on  a  log  scale  with  time  as  a  parameter. 

Graphical  solution  of  time-temperature  history  of 
bakelite  surrounding  a  copper  plug  while  the  plug 
is  being  cooled  by  the  high  velocity  gas  stream. 


13  Temperature  above  adiabatic  wall  versus  time , 

Comparison  of  graphically  predicted,  and  experimentally 
measured  temperature -time  histories  for  two  points  in 
bakelite  surrounding  a  copper  plug. 
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FIGURE  1 1 


ABSTRACT 


The  objectives  of  this  program  are  to  de,termine  the  local 
convection  heat  transfer  coefficients  in  the  convergent  portion 
and  throat  sections  of  de  Laval  nozzles;  and  to  determine  the  ef¬ 
fect  on  the  local  heat  transfer  coefficients.,  part'icularly  at  the 
throat,  of  varying  the  contour  in  the  converging ‘section  of  the 
nozzle. 

Three  two-dimensional  nozzles  (described  hereafter  as  III, 

IV,  V),  with  pressure  gradients  as  defined  in  Fig',  la,  and  result¬ 
ing  cross-sectional  areas  as  described  in  Fig,. lb,  have  been  tested. 
Heat  transfer  coefficients  have  been  obtained  for  both  the  flat  and 
contoured  walls. 

The  results  do  not  correlate  when  using  a  length  dimension 
based  on  hydraulic  diameter,  as  conventionally  employed  for  estab¬ 
lished  pipe  flow. 

The  results,  for  both  the  flat  and  contoured  walls,  do  corre¬ 
late,  when  using  a  length  dimension  measured  from  the  nozzle  entrance 
along  the  walls,  .however,  the  heat  transfer  coefficients  are  higher 
than  for  turbulent  boundary  layei  flat -plate  flow.  This  appears 
to  be  partially  due  to  the  technique  used  in  obtaining  the  appar¬ 
ently  local  coefficients ,  and  par  tially  due  to  the  negative  pressure 
gradient  effect, 

A  significant  result  of  these  tests  is  that  the  throat  coeffi¬ 
cients,  for  the  same  mass  flow  rates  through  all  the  nozzles  III, 

IV,  and  V,  agree  within  t  15^,  indicating  that  the  markedly  dif¬ 
ferent  throat  approach  contour  variations  shown  in  Fig,  2  do  not 
produce  significant  variations  in  the  nozzle  throat  coefficients, 

A  theoretical  solution,  predicting  the  local  heat  transfer  co¬ 
efficients  along  the  flat  wall  of  nozzle  III,  has  been  made  ^y 
troducing  the  pressure  gradient  of  nozzle  III  into  the  von  Karman 
boundary  layer  momentum  equation.  The  results  of  this  work  are  com¬ 
pared  to  the  experimental  findings,  and  are  utilized  to  interpret 
the  experimental  results , 
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PREDICTION  OF  THE  TIME-<CEMPERATURE  HISTOir  CF  THE  THERMAL 
PLUG,  WHILE  COOLING,  USING  A  GRAPHICAL  TECHNIQUE 


In  the  process  of  attempting  to  analyze  the  time -temperature 
history  of  the  thermal  plug  correctly,  an  approximate  method  was 
derived  for  evaluating  the  time -temperature  history  of  a  body 
which  losses  thermal  energy  via  two  paths  such  as  the  plugs  in 
this  experiment. 

The  method  will  be  described  by  using  the  plugs  as  a  basis 
of  analysis  but  it  can  readily  be  seen  that  the  method  is  adap¬ 
table  to  any  similar  system. 

By  making  the  following  assumptions,  the  cooling  time- 
temperature  history  of  a  plug  can  be  predicted  as  a  function 
of  plug  dimensions  and  the  gas  film  heat  transfer  coefficient. 

(a)  Assume  time  of  heating  plug  and  temperature  to  which 
it  is  raised  are  known.  Then  the  temperature  distri¬ 
bution  in  the  bakelite  at  time  zero  of  the  cooling  run 
is  known. 

(b)  Assume  radial  heat  flow  only  from  the  plug  into  the 
bakelite . 

(c)  Assume  the  thermal  dlffusivity  of  copper  with  respect 
to  heat  flow  from  its  face  to  the  gas  stream  and  from 
its  side  area  to  the  bakelite  is  sufficiently  high 

to  consider  the  plug  at  isothermal  conditions . 

(This  has  been  verified  for  the  case  of  the  experimental 
data .) 

(d)  Assume  the  contact  resistance  between  the  copper  and 
bakelite  is  negligible. 

(©)  Assume  heat  flows  from  the  plug  only  via  the  two  paths 
just  mentioned. 
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if  r  is 


Now  can  be  v/ritten  as  t^,  - 

b  A  r 

small  as  compared  to  r  which  is  the  case  here. 

First  consider  the  time-temperature  cooling  history  of  the 
plug  when  the  stream  face  is  insulated,  =  0, 

Then  from  a  consideration  of  rate  equation  and  energy  balance 
on  the  plug  for  a  very  short  time  in  finite  differential  form 


V/DC  Atg  = 


Va. 

4r 


(tg  -  tj,)as 


A  r  tg  - 

^b^  A  ®  ^ 


Consider  the  graphical  solution  described  earlier,  of  the 
time- temperature  history  of  the  bakelite  while  the  plug  is  heating. 
It  is  solved  for  the  pl.ug  heating  time  and  a  sample  solution  ap¬ 
pears  in  Fig,  (11),  If  some  value  A  a  is  chosen  such  that 


~  ^r 

A  a  "  A  r 


then  equation  (2)  results. 


~  ^r  ^  Al 
A  tg  A  a 


(2) 


Substituting  equation  (2)  ^n  equation  (1) 


V^C 

V2 


Ar 

A9 


=  AL 
A  a 


k,  A.A0 
Aa  =  -AA 


V^C 

This  is  assuming  is  constant  for  time  A  0. 

Actually,  if  A  a  is  small  compared  to  A  r  this  is  very  nearly  true 
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Since  all  the  physical  properties  other  than  temperature 
are  previously  related  or  known  in  order  to  make  the  graphical 
solution  of  the  time-temperature  history  of  whe  hakelite  while 
the  plug  is  heating*  it  remains  to  introduce  ^  a  into  the  graph¬ 
ical  solution  to  predict  the  time-temperature  history  of  the 
plug  during  the  cooling  run.  By  laying  it  out  as  shown  in  the 
following  figure  tg'  which  equals  ^  ^ 

0  0 

directly  and  the  time-temperature  history  of  the  plug  and  bake 11 te 
matrix  can  be  predicted  in  rapid  fashiono 


In  the  case  where  there  is  heat  loss  to  the  bakelite  and 
gas  stream  simultaneoasly,  divide  the  mass  of  the  plug  into  two 
components  and  Y/o)  Let  one  supply  thermal  energy 

to  the  gas  stream  and  one  supply 

thermal  energy  to  the  bakelite » 


Then  from  consideration  of  thermal  energy  flow  rates  and 
energy  balance  on  the  plug  the  following  two  equations  are 
obtained 8 


V/O  )3,  c 


A  t 


Ml  ^^0  "  ■‘^aw^ 


A  t 


0 


Ao 


tg  -  t 
=  k^A  (— ^ - 


(5) 


The  A©  to  be  used  is  evaluated  from  the  physical  proper 
ties  and  the  selected  value  of  A  r  for  the  graphical  solutiono 
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A  tg  must  be  the  same  in  equations  (4)  and  (5).  Then  solving 
for  A tg  in  both  equations  and  equating  gives 

A  C 


'*0  -  ^aw>  = 


^b^2 


(^0  - 
Ar 


(6) 


or  -  ^aw  .  II.  Vg 

^0  -  ^  ^2 


(7) 


ka. 


h  P 

V,  k  A  -n  ~  if  which  is  a  constant  for  any  cooling  run, 


Then  equation  (7)  becomes 


~  -V-  =  ( f-  - 
^0  ■*  >  ^2 


where  ^  +  Vg 


or 


^2  = 


1+1- 

(T 


t 


9  aW' 
*0  -  ’^r' 


(7a) 


Substituting  equation  (5)  in  equation  (7)  gives S 


At 
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.  Vg^s 


.  Vg^  ® 


*0  ■  *r  1^)  C<ar  /°C  VAr 
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,  ,  1  ^^0  "  ‘‘^aw^l 

1  +  ~  Tz - r-v“J 


Urn  W  W*  ff 

3'  (tg  -  t^) 


(8) 


Letting  A  t^  _  '^9  "  "^r 
A  a  A 

A  r 
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(9) 


or  A  t, 


Substituting  equation  (9)  in  equation  (8)  gives 


Aa  _  ^ ^ 

“  V/OC 


^^0  ' 

-  tj.) 


]  - 


(10) 


Letting 

yAav  =  constant 


^  (t!  ■.  t  T  J 


(11) 


Now^  ,  ^5  t^^,  r  are  all  constants o  (tg  -*  t^)  and 

(tg  “  can  be  read  from  the  graphical  solutiono 

To  predict  time -temperature  cooling  history  of  a  plug 
then  the  procedure  would  be  exactly  the  same  as  in  the  case 
where  the  gas  stream  face  of  the  plug  was  insulated  except  that 
in  this  case  /t „  is  not  constant  but  will  have  to  be 


evaluated  each  time  as  a  function  of  tr,  -  t 

6  aw 


It  is  seen  that  A  a 

Ar 


t^  “  t 
9  aw 


will  vary  as  a  function  of 


In  this  case  (tg  -  t^^^)  will  decrease  less 


rapidly  than  (tg  -  t^) 


Therefore  A  a  will  approach  a 
Ar 


very  large  quantity  and  the  solution  will  be  invalid  in  this 
reglono  Hov/ever,  as  -  t  )  approaches  zero  the  heat 


•»  3l  - 


loss  to  the  bakelite  becomes  negligible  and  ^  tg  can  be 
evaluated  in  this  region  by 


Ato  hnA, 

_ 0  _  1  X  /  i 

A9  "  V/oC 


) 


This  method  was  used  to  repredict  the  cooling  time- 
temperature  history  using  the  heat  transfer  coefficient  de¬ 
termined  by  the  method  described  earlier  in  the  appendix. 

The  results  are  shown  in  Fig,  (l4).  The  repredicted 
time-temperature  history  shows  very  good  comparison  with  the 
measured  time-temperature  history  of  the  plug  and  it  is  con¬ 
cluded  from  these  results  that  both  the  method  just  described 
and  the  method  of  determining  the  heat  transfer  coefficient 
are  valid. 
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Figure  iH-  Temperature  of  copper  plug  above  adiabatic  wall 
temperature  versus  timeo 

Comparison  of  the  experimentally  measured  time- 
temperature  history  of  the  plug  and  the  time- 
temperature  history  predicted  by  a  graphical 
technique o 
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TEMPERATURE  OF  PLUG  (Tp- 


PLUG  COOLING  TIME  SECONDS  FIGURE  1  4 


INVESTIGATION  OF  THE  THERMAL  PLUG  EFFECT 
IN  ESTABLISHED,  TURBULENT  PIPE  FLOW 


To  obtain  more  quantitative  experimental  evidence  of  the 
significance  of  the  delayed  thermal  boundary  layer  with  respect 
to  the  hydrodynamic  boundary  layer  a  plug  identical  in  geometry 
to  those  used  in  the  nozzle  experiments  was  placed  in  the  wall 
of  a  pipe,  insulated  from  the  pipe  by  a  bakelite  annulus,  and 
both  the  bakelite  and  plug  were  machined  flush  with  the  inner 
surface  of  the  pipeo  The  same  transient  technique  was  used  to 
obtain  the  heat  transfer  coefficient  in  this  test  apparatus  as 
was  used  in  the  nozzle  experiments  <> 

The  flow  rate  was  sufficiently  high,  Re^  =  2o44  x  10  ,  and 

the  pipe  length  great  enough,  L/D  of  70,  to  assure  established 
turbulent  pipe  flow. 

The  predicted  Stanton  Noo  based  on  the  correlation 
St  Pr^/^  =  o023  Rep" 
was  Sto  NOe  =  2o46  x  lO"^  o 

The  exuerimental  St,  No,  value  obtained  using  the  thermal 
plug  was  o 

St,  No,  =  5,70  X  10“-^ 

If  Stanton’s  No,  were  predicted  on  a  flat  plate  basis 
using  the  maximum  velocity  of  the  pipe  as  the  free  stream 
velocity  for  flow  over  a  flat  plate  and  the  length  in  the 
Reynolds  No,  as  1/2  the  length  of  a  square  plug  equivalent 
in  area  to  the  round  plug  tested,  the  Stanton  No,  value  based 
on  the  correlation 

St  Pr^/^  =  o0295  Re^  " 

Jv 

would  be  St,  No,  =  6,25  x  lO"  , 

It  was  anticipated  that  the  experimental  results  would 
lie  between  the  Stanton  No,  predicted  for  established  flow 
and  established  thermal  boundary  layer  and  the  Stanton  No, 
predicted  on  a  flat  plate  basis  considering  both  the  thermal 
and  hydrodynamic  boundary  layers  as  starting  at  the  leading 
edge  of  the  thermal  plug. 
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The  conclusions  that  can  he  drawn  from  the  experimental 
results  are: 

1)  The  anticipation  of  the  results  lying  between  the 
two  limits  just  described  is  valid, 

2)  That  the  results  are  nearer  the  upper  limit 
St  =  6  .25  X  10 . 
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